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1 Introduction 



The nonlinear Schrodinger (NLS) equation is one of the universal nonlinear integrable equation since 
it describes the slow modulation of the amplitude of a weakly nonlinear wave packet in a moving 
medium pQ. It has been used with great success in nonlinear optics to describe the propagation 
of sufficiently broad pulses, or slowly varying wave trains whose spectra are narrowly localized 
around the carrier frequency [2J. Today high-speed fiber-optic communication demands ultra-short 
pulses and technological progress for creating them was achieved, however, in these conditions the 
description of the evolution of these pulses lies beyond the usual approximations leading to the NLS 
equation. In [3j Schafer and Wayne proposed an alternative model to approximate the evolution 
of ultra-short intense infrared pulses in silica optics. After a scale transformation of variables their 
short pulse (SP) equation can be written as 

Uxt = u + ~ {u 3 ) xx or u t = (d^u) + ^u 2 u x , (1) 

with u = u(x, t) representing the magnitude of the electric field and subscripts x and t standing 
for partial differentiations. It was shown in [2] by numerical simulations that the SP equation is a 
success to describe pulses with broad spectrum. 

In the last years the SP equation became a comprehensively studied equation of soliton theory. 
The integrability of this equation was first established in [5] where for the linear spectral problem 
of the Wadati-Konno-Ichikawa type 

<5> x = X$> , $ t = T$ , (2) 
the corresponding zero-curvature representation (ZCR) 

X t -T x + [X,T] = 0, (3) 



is given by 



A Xu x \ ^ / fu 2 + 3 L | {u s ) x - \u ^ 



\u x -X J \ I (u 3 ) x + \ 



_A U 2 



2 " 4A 



Also in [5] a transformation was found which relates the SP to the sine-Gordon equation and it 
was used in [6] to derive exact soliton solutions of the SP equation from known soliton solutions 
of the sine-Gordon equation. In [7] the SP equation was studied as one of four Rabelo equations 
of differential geometry [8]. The recursion operator, Hamiltonian structures and conservation laws 
for the SP equation were also found in [9l [10]. The first Hamiltonian structure follows from the 
Lagrangian density for the SP equation in nonlocal form (prj 
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and reads 

Aff„ r r i i . _i 

(5) 



8H2 

ut = T>\ — — with T>i = d , H2 = I dx 

ou 



1 4 -*-/«— 1 \2 

— 7 U [O It) 

24 2 V ; 



The second Hamiltonian structure follows, for instance, from the recursion operator obtained in [5] 
using the cyclic basis technique [TTJ [12] and yields 



T>2-^- with T>2 = d 1 + u x d 1 u x , Hi = - [ dxu 2 . (6) 
ou 2 J 



Many other results were found concerning the SP equation such as Hirota's bilinear representation 
J, multisoliton solutions [15J and periodic solutions [HUE], among others. 



To take into account effects of polarization and anisotropy (see [18] and references therein) 
two-component (u, v) integrable generalizations of the SP equation were proposed in the literature 
[19|-|23j. These equations play the same role in dynamics of ultra-short pulses as the pair of 
coupled NLS equation does for broader pulses [25]. Note that some of these systems reduce to 
the SP equation ([T])ifii = 0or7; = while other systems if u = v. The integrability of these 
systems were obtained mainly from a zero curvature representation or bilinear formalism. An 
algebraic approach to the integrability based on a Hamiltonian formulation for infinite-dimensional 
dynamical systems [26] is lacking and in what follows we address the bi-Hamiltonian integrability 
of some of these systems written in evolutionary form 

\vj t \6H 2 /Sv J ySHt/Sv J 

As is well known from the Hamiltonian structures T>\ and T>2 we can construct a recursion operator 
R = T>2T>^ 1 that yields the hierarchy of equations, symmetries and conserved charges associated 
with these two-component SP equations [27j . 



2 Pietrzyk— Kanattsikov— Bandelow Systems 

Pietrzyk et al. [19] introduced the following three integrable two-component SP equations 

u xt = u+- (n 3 + 3uv 2 ) xx , v xt = v + - (v 3 + 3u 2 v) xx , (7) 

u xt = u + - (u 3 - 3uv 2 ) xx , v xt = v - ^ (v 3 - 3u 2 v) xx , (8) 

I 1 

u xt = u + - (u 3 ) xx , v xt = v + - (u 2 v) xx , (9) 
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by a direct generalization of the SP equation's ZCR matrices They follow from ([3|) with 4x4 
matrices given by 



X 



AIL, -AI 
where I is the 2x2 identity matrix and 



T 



§u 2 + ii |(u 3 ) x -iu 
|(u 3 ) x + iu 



-|u 2 



(10) 



u 



?; it 



u 



u v 



u 



reproduce ((7J), (|HJ) and ([9]), respectively. In [2l] Sakovich pointed out that equations (J7J and (jHJ) 
can be reduced to two decoupled SP equations ([T]) in the new variables (p = u + v , q = u — v) 
and (p = u + iv,q = u — iv), respectively. So, ([7]) and © are not a short pulse counterpart 
of the Manakov system [25J of coupled NLS equations where the polarization modes do interact 
nonlinearly |28| . For completeness we will obtain the bi-Hamiltonian structure for (J7J and ([8]). 
These equations written simultaneously in a nonlocal form 



ut = d l u + 



6 



u 



± 3ra 2 ) 



v t = d~ l v ± - (v 3 ± 3u 2 v) 



follow from the following Lagrangian density 



£ = i; ((TV) - \u 3 v t ^v 3 n + {d~ l u) (d~ l v) , 



6 6 

and we can read the first Hamiltonian structure directly, 



d 
d 



H 2 



dx 



-u 3 v ±-v 3 u - (d- l u)(d- l v) 
6 6 



(12) 



From the diagonal second Hamiltonian structure for the decoupled system with elements given by 
([U]) in the variables p and q 



d 1 +p x d 1 p x 

cT 1 + q x d^q x 



1 



dx (p 2 + q 2 ) . 



we use the usual results [29J for transformation laws of Hamiltonian structures under the above 
transformations (p, q) — > (u, v) to get 



d 1 + u x d l u x ± f x 9 1 v x u x d l v x + v x 5 u a 



u x d 1 v x + f x 9 1 u x 



#1 



dx (u 2 ± v 2 ) . 
(13) 
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The proof of the Jacobi identity for these structures (and the next ones) as well their compatibility 
necessary to satisfy Magri's theorem [30] can be proved by the standard method of prolongation 
described in |27] in a straightforward way. 



Equations ([9]) written in a nonlocal form 



Vt 



d 1 v + - (u 2 v) x , 



(14) 



can be derived from the following Lagrangian density 



1 



C = v (d- l ut) - -u 3 v + (d^u) (d^v) , 



which gives the first Hamiltonian structure 

Mao)' "> 



dx 



-u 3 v - (d^uYd^v) 
6 



(15) 



Also, was remarked in [I]3] that the system §§§ describes the propagation of a small perturbation v 
on the background of a solution u of the scalar SP equation (fTJ) and from this perspective we can 
rewrite Q as 



u x t = u + - {u 3 ) xx = K(u) , v xt = K'(u)[v] 



(16) 



where 



K'(u)[v] = —K{u + ev) 



e=0 



(17) 



is the Frechet or directional derivative of K(u) in the direction of v. So, the equation for v is the 
linearized equation for u given by ([T]) and objects for the system ([9]) or (fl~6j) such as conserved 
charges and Hamiltonian structures can be obtained directly from the corresponding objects of the 
scalar SP equation (fTJ) using the formulas given in [31] such as 



V 



V 
V V'(u)\v] 



H = H'(u)[v] 



where the "~" referes to the objects of the scalar SP system ([T]). Hence, the second Hamiltonian 
structure for (|14p can be obtained from ([6]) (as well the first one f)15|) from ([5])) and it is 







d + u T d 



li T v. 



d + u x d l u x 
: d~ 1 u x + u x d~ 1 v a 



/ 



dx uv . 



(18) 
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3 Dimakis— Mmler-Hoissen— Matsuno System 



From a bidifferential approach to the AKNS hierarchies [20] or from a bilinear formalism combined 
with a hodograph transformation [21] these authors proposed the multi-component system 



1/ ^ A 

Ui,xt = Ui + ~ \Ui tX u j ] ! 1 = 1,2, ...n 



2 . 

with reduced two-component SP system given by 

u xt = u+^ [(u 2 + v 2 )u x ] x , v xt = v+ - [(u 2 + v 2 )v x ] x , 

where u = u\ and v = v>2- After the transformations (u,v) —> ((u + v)/2, (u — v)/2i) this system 
can be written as 

u xt = u + ^ (uvu x ) x , v xt = v + - (uvv x ) x , (19) 
which is integrable with zero curvature ([3]) with matrices given by 

uvu x — \u \ 

(20) 





f A \u x \ 




H 








\ \v x -A J 





^UVV X + \v — IjUV — 

A Lagrangian density for the equations (|19p written in a nonlocal form can be easily obtained and 
it is 

C = uv t - uid^v) - \u 2 (v 2 ) x , 

o 

and we can read the first Hamiltonian structure as 

(21) 



P, - ( J ^ ! . i/ 2 / d.r 



u{d- 1 v) + \u 2 {v 2 ) x 



To find the recursion operator of the system (|19p . we use the cyclic basis method [11} I12j, We 
pose the problem of finding all the evolutionary systems 

ut = f, vt = g, (22) 

possessing the ZCR ([3]) with the matrix X given by (|20p and the matrix T being initially not fixed. 
We rewrite the ZCR in its equivalent (covariant, or characteristic) form 

fC tt + gC v = V x S, (23) 
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where C u and C v are the characteristic matrices, 

'ax' 



C u — ~ V a 



dUn 



C v — — V 3 



ax 

dv x 



(24) 



since X = X(u x ,v x ). The covariant derivative operator V x is defined as X? X H = dH — \X,H\ for 
any 2x2 matrix H, and the matrix S is related to the unknown matrix T as 



Q T f 9X dX 



(25) 



du x dv x 

For X given by (|20p we find that the cyclic basis is three-dimensional and consists of the matrices 



C u 



-\ 2 v x 2A 2 



AV 



C r 



X 2 u x 
-2A 2 -X 2 u Xi 



2\ 3 (2 + u x v x 
2\ 3 v x + \ 2 v x 



2X 3 v 2 x 

and the closure equations for the basis being 

V X C V = a\C u + a2C v + a^V X C U , 
V 2 X C U = hC u + b 2 C v + b 3 V x C u , 

where 



ai = A 
bi = -A 
b 2 = -A 
= -A 



(2 + u x v x )u x 



V X U XX ~\~ 2u x v xx 



, a 2 = A 
2 -I - u x v x )v xxx 



(2 + u^J 



2 

V X U XX 



+ X< 



u 



03 



(2 + u x v x )v 2 u xx 2 



3v x v x 



+ X 2 



V -^L - (2 + u x v x )v 2 x 



(2 + u x v x ) 



+ 



We decompose the matrix S over the cyclic b 

S = pC u + qC v + r V X C U , 



(26) 



(27) 



(28) 



(29) 
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and we obtain from ()23[) and (|27p that 



/ = dp + a\q + bir , 
g = dq + aiq + b 2 r , 

P = -9r - a 3 q - b 3 r , (30) 

which can be rewritten, using (|28p . as 

Q = (M + AL + A 2 ^) f j , (31) 

where the 2x2 matrix differential operators M, L and K do not contain the spectral parameter 
A, and their no null components are 

Mn = —9 + _ ; Mi 2 = -a 2 - ^££3 + ^ - . Af 2 i = d , 

V xx V X x V xx V X x V xx V xx 

t 2 ^X^xj^xx V X U XX V X V XXX 

L n = u x , L 2i = 2 + u x ?j x , L22 = ~3v x v xx H , 

Uxx ^ XX ^ XX 

T f n , V X U XX \ _ ^"rai . / n V X U XX \ f xxx 

L12 = 12 + UajWj. I d - 2t> x u xx - h I 2 + n x -u x H I , 

\ "ra / V xx \ V xx J V xx 

K 12 = -ulvl + i;1 ' , i^ 22 = -(2 + u x v x )vl + r ^—EL . (32) 

The hierarchy of evolutionary systems (|22p possessing ZCR ([3]) with the matrix X given by (|20p 
is completely determined by the relation (|3ip . where q and r must satisfy the conditions df /dX = 
dg/dX = 0. If we would have K = in (|3ip . we could immediately write down the recursion 
operator for the represented hierarchy as R = ML~ l [12] . The more complicated case of the 
relation (|3ip with .fT 7^ was studied in [13], where it was shown that the recursion operator R is 
given by 

R = MN' 1 , N = ML~~ X K - L , (33) 
if the operators L and K satisfy the condition 

KL~ l K = . (34) 
In the present case of X it turns out that the operators L and K do satisfy the condition (|34p . and 
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we have from (|32p and (|33p the following: 



/ a —a — av T 



M 



^ (2 + ti^)— - u 2 -(2 + u x v x )—dv xx + m^j^ 



N 



-(2 + u^) 



+ 



-SUaja, + V x V a 



(35) 



and then 



|a - \dFu x d~ l Fv x d -\dFu x d- x Fu x d 



R 



\dFv x d~ 1 Fv x d 



\d + \dFv x d~ l Fu x d 



where 



f = (i + ^r 1/2 . 

Its inverse can be easily obtained to yield 

ld~ x + u x d~ 1 v x —u x d~ 1 u x 

v x d~ 1 v x —2d~ 1 — v x d u x 



R- 



and from the factorization R 1 = T>2'D 1 1 and the first structure (j21h we get the second Hamiltonian 
structure of the Matsuno system 



u x d 1 u x 2d 1 + u x d 1 v x 
2d~ 1 + v x d~ 1 u x v x d~ 1 v x 



H i = - / dxuv . 
2 



4 Feng System 



Using two sets of bilinear equations of a two-dimensional Toda lattice linked by a Backlund trans- 
formation 

r D s D y f ■ f = ±(/ 2 - f) , r £>,£> y s • 5 = I( 5 2 - g 2 ) , 
D S ZV • / = ±(/ 2 - / 2 ) , (D s D y g .g=l(f- g 2), 

a particular hodograph transformation 

(x = y-(ln(FF)) s , 
It = s , 
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and the dependent variable transformation 

u = i ( In ^ , v = i ^ln ^ 

where i 7 = fg, G = fg and F and G are the complex conjugate of F and G respectively, Feng |22] 
proposed the coupled SP equation 

Wxt = u + uuj + i (u 2 + t> 2 ) , v xt = v + vvl + - (u 2 + u 2 ) . (36) 

From a prolongation study [32J of this system we have obtained for the zero curvature ([3D the 
matrices 



X 



A(l + u x v x ) X(u x - v x ) 
\{u x -v x ) -\{l + u x v x ) 

|(u 2 + v 2 ){l + u x v x ) + i |(u 2 + v 2 )^ - Ua .) - 1(« - v) 



(37) 



| (u 2 + v 2 )K - « x ) + i(u - u) — |(u 2 + ^ 2 )(i + «x«x) - IX 

Equations (|36p written in a nonlocal form can be derived from the following Lagrangian density 

Is Is 

£ = u xt v — uv v u xx u v xx , 

6 6 

and we can read the following first Hamiltonian structure 

£>i = ( g° i ^ ^ , H<2 = j dx ( uv + \ y3 ' Vxx + \ v * u **\ • ( 38 ) 

It is not possible to obtain a recursion operator from the matrix X in (|37|) by cyclic basis 
technique, for the same reason as the Example 6 in Namely, the transformation (|43p - (j45p . 

given below, changes the linear problem ([2]) with the matrices (|37p into the linear problem ty y = Y^/ 
and \£ s = S*&, where *$>(y,s) = <3?(x,i), with the matrices 

Y _ fXcosa' Asino-' \ _ / j\ ~2 a s\ 

' "VAsina' -Acos^' [±a' s ) ' 

where two dependent variables u and v have merged into one dependent variable a'. However, 
we can derive the Hamiltonian structures for the Feng system from the well-known ones of the 
sine-Gordon (SG) equation a ys = sin a, in light cone coordinates [33] , 

P (ff) = dy 1 , F (ct) = - fdy cos a, (39) 
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and 



= 8 y + (Jyd^cjy, = ~ J dy cos (J ( d~ 1 sin <r) 2 , (40) 

using the transformation (x,t, u, v) — > (y,s,a,a r ) between the Feng system and the pair of uncou- 
pled SG equations, a ys = sincr and a' ys = sincr', discovered in [22]. In fact, the Feng system (j36j) 
and the pair of SG equations in evolutionary form 



where 

/ = 2 + ' L ' 2 ) u * + ®x Xt L u ~ vu x v x ), 
g =^(u 2 + v 2 ) v x + d' 1 ^ - uu x v x ), 
h = dy 1 sincr, /i' = dy 1 sin a' , 
are related to each other by the Feng's transformation 



a = arccos — h arccos -, (42) 

p q 
l l 

a = arccos arccos -, (43) 

V Q 

t = s , x = z(y,s) : z y = —, (44) 

pq 

p 2 = l + u 2 x , q 2 = l + v 2 x . (45) 



Under the transformation (|42p ~ (|45p the right-hand sides of evolution systems (|4ip are related 
to each other as 

with a 2 x 2 matrix linear operator P. Applying d s to (|42p and (|43p . we get 

1 1 

h = -^{u xx z s + d x f) + -?{v xx z s + d x g), 

(47) 

11 

/i' = -f(Ux X Z s + C^/) Av xx Z s + <9 X .G/). 
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Applying d s to the relation z y = l/(pq) of (|3S|1 and taking into account that 

1 



we get 



d v = ~ 9 ^ d y 1 = d x 1 pq, 



*s = ~d x 1 [K + u x )/i + («x - v x )ti] . (48) 



Eliminating z s from (|47p and (|48p . we obtain (|46p with 



i /V + u X x3 x 1 (u :c + Wa;) p 2 + u xx d x 1 (u x - v x ) 
P=^d- 1 l | . (49) 

\ Q 2 + v xx d x 1 (u x + v x ) -q 2 + u^S,. 1 (u x -v x ) 



2 



Let us consider the systems (I4ip as Hamiltonian ones, that is 

Jj _ qj(u,v) ^(u,v) jj(u,v) _ jy(<r,<r') ^(cr.cr') jj{<j,<j') 

where 



P are Hamiltonian operators, H are Hamiltonian functionals, and their superscripts (u, v) and 
(a, a') indicate which system they correspond to. We have already found that Ut = -PS S with P 
given by (|4"9"]) . Let 



= PV^Q, (50) 

so that 

A (u,v) H (u,v) = Q-l A (*,<r') H (a,a>)^ ^ 

with a 2 x 2 matrix linear operator Q. From the relations 

= d t H^= [dxUTA^H^ = fdy — (P-E s ) T A^H^= fdyT^—P^A^H^ u ' v \ 

J J pq J pq 

= d s H^ a ^ = jdyT^A^'^H^ = J dyY^QA^ H^ v \ 
where the superscripts T and f stand for the matrix transposition and operator conjugation, re- 
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spectively, we see that we should set 

Q = 

pq 

Due to f4*9|) . we have 

-1 /V - ( n 2 + v x )d~ 1 u xx q 2 - (u x + v x )d~ l v x 
2pq \P 2 ~ (u x - v x )d~ l u xx -q 2 - (u x - v x )d- l v. 
and the explicit expression 

/ 



Q 



d x \ 



-i 



Q 



will also be useful. 



-d x 



P 

p 
\1 



1 + ( + -J 



1 + v x d„ 1 + 



1 
P 

P 

q 



1 + u x d„ M — 



1/ ^ieie "xx 



-l + u-a. 



1/ l^xx V xx 



(52) 



(53) 



Now, using the relations (1501) and (ISTT) with (|49|) . (I52j) . and (j53j) . we can derive the Feng system's 
Hamiltonian structures from the ones of the scalar SG equation. We can take 



2?W N 
j 92?( ff/ ) 



ff foe') = ijy(«r) + L H (v') 



a 



(54) 



for the pair of uncoupled SG equations, where a and /3 are arbitrary nonzero constants, while T> 
and H with the superscripts (a) and (a') are given either by (f39l) or by fTOl) . Choosing ([39]) . we 
have 



I (3d: 



For the Hamiltonian functional, we obtain the following: 



— sm cr 



i sin <t' 



/J 



(55) 



I _ 1 ) ( ^ 

/3 ay 



1 1 



dx u T v T . 
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The choice of a = f3, which would seem most natural, is actually not allowed and without loss of 
generality we set 

a = _2, = 2 (56) 

to get 

H x = #(«•") = - J dx u x v x . (57) 

For the Hamiltonian operator, we obtain from ([55]) with ([56]) via ([50]) . ([4*9]) . and ([52]) the following 
expression: 

',d x x b — bd x x a) b d~ x d — a d~ x 
V 2 = P^'") = I I , (58) 



where 



a = d x 1 u xx d x 1 v x , b = d x l + u x d x 1 u x , 
c = d~ x v xx d~ x u x , d = d x l + v x d x l v x . 

(bv xx ) = (v x + c^u xx ) , (du xx ) = (u x + a^^) , 



Using 



it is straightforward to check that (j36|) in evolutionary form follows from (|57[) and (j58p as a Hamil- 
tonian system. In completely the same way, using the SG equation's Hamiltonian structure (j4"0l 
and the choice (1561) in (1541. we obtain (1381). 



5 Yao— Zeng System 

Yao and Zeng [23] claim to have constructed a new spectral problem ([2]), with matrices X and T 
spanned by a specific loop algebra, which generates a new hierarchy of coupled SP equations having 

Uxt = u + - {u 3 ) xx , v xt = v + ^ (u 2 v x ) x , (59) 

as its first member. However, the ZCR obtained in their paper does not reproduce (|59j) as can 
be easily checked. Also, by a direct analysis of their ZCR, one can prove that the system (|59p 
cannot be associated with their spectral problem. Here we establish the integrability of (|59[) in the 
following way. By the transformation v — > v x the Pietrzyk et al system ([9]) can be mapped into 



14 



(|59p and particularly the ZCR matrices (jlOp yields, after a gauge transformation [34J, 



X 



( 


A 


Alia; 





o \ 




Xu x 


-A 













Auje 


A 


Xu x 


V 







Xu x 


-A / 



/ -4 + 



r 



A. ,2 



U U x + t;U 



4A 



\ — -|n 2 t; x — 












4A + 2 u2yjx 2 U 



This corrects the ZCR and in terms of the Lie algebra spanned by the generators {ei, e 2 , e3, e4, es, ee} 
considered in [23] we have the following decomposition 



X = Aei + An :c (e2 + e 3 ) + Xv x (e 5 - e 6 ) , 
T = Tnei + Ti 2 e 2 + T 2 ie3 + T 32 e5 + T^ee , 

in disagreement with the one provided by Yao and Zeng in their isospectral problem. 



For completeness, the bi-Hamiltonian structure of the Yao-Zeng system (|59p in evolutionary 
form can be obtained from (115p and (|18p using the transformation laws of Hamiltonian structures 
under the map (u,v) — > (u,v x ), and they are 



-1 

1 



H 2 



dx 



1 ., , 
-u v x — (d u)v 



6 



-d 2 - -u^d 1 

<9~ 2 + d~ 1 u x d~ 1 u x —d^ 1 {v xx d~ l u x — Uxd^Vxx) d^ 1 



I 



H, = dx UV r ■ 



In [3l] we provide details of how the Yao-Zeng system (|59|) is related with a perturbation of 
the SP equation ([!]). 
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6 Conclusions 



In this paper, we have obtained the bi-Hamiltonian structure for some of the two-component SP 
equations proposed in the literature to take into account the propagation of polarized pulses in 
anisotropic media. Namely, we studied systematically the Pietrzyk-Kanattsikov-Bandelow [19j, 
Dimakis-Miiller-Hoissen-Matsuno [20\ I21j. Feng |22j and the Yao and Zeng [23] systems. Also, we 
have obtained the zero-curvature representation for the Feng and Yao-Zeng equations. 

Through a Painleve analysis of two one-parameter equations belonging to the general two- 
component short pulse equations for cubically nonlinear anisotropic optical fibers, introduced in 

HSI, 

Un,xt — CniUi ~\~ C n ijk{UiU ~jUk)xx j i>i3i ^ — 1 5 2 , 

where c n i and c^.-^ are constant coefficients determined by optical properties of the medium, two 
more integrable systems besides ([7])-([9]) were found in [24] 

u xt = u+- (u 3 + uv 2 ) xx , v xt = v + - (v 3 + u 2 v) xx , (60) 

U °* = U + 6 ^xx > V *t = V + 6 ( U2y )xx . ( 61 ) 

where u and v denote the polarization components U\ and U2- No ZCR and bi-Hamiltonian 
structures are known at this point and these systems are still under investigation. 
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